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ABSTRACT 

We investigate the radial distribution of galaxies within their host dark matter halos as measured in 
the Sloan Digital Sky Survey by modeling their small-scale clustering. Specifically, we model the Jiang 
et al. (2011) measurements of the galaxy two-point correlation function down to very small projected 
separations (10 < r < 400/i~ 1 kpc), in a wide range of luminosity threshold samples (absolute r-band 
magnitudes of —18 up to —23). We use a halo occupation distribution (HOD) framework with free 
parameters that specify both the number and spatial distribution of galaxies within their host dark 
matter halos. We assume one galaxy resides in the halo center and additional galaxies are considered 
satellites that follow a radial density profile similar to the dark matter Navarro-Frenk- White (NFW) 
profile, except that the concentration and inner slope are allowed to vary. We find that in low 
luminosity samples (M r < —19.5 and lower), satellite galaxies have radial profiles that are consistent 
with NFW. M r < —20 and brighter satellite galaxies have radial profiles with significantly steeper 
inner slopes than NFW (we find inner logarithmic slopes ranging from —1.6 to —2.1, as opposed to —1 
for NFW). We define a useful metric of concentration, M 1 / w , which is the fraction of satellite galaxies 
(or mass) that are enclosed within one tenth of the virial radius of a halo. We find that Mx/xo for 
low luminosity satellite galaxies agrees with NFW, whereas for luminous galaxies it is 2.5 — 4 times 
higher, demonstrating that these galaxies are substantially more centrally concentrated within their 
dark matter halos than the dark matter itself. Our results therefore suggest that the processes that 
govern the spatial distribution of galaxies, once they have merged into larger halos, must be luminosity 
dependent, such that luminous galaxies become poor tracers of the underlying dark matter. 
Subject headings: cosmology: theory — galaxies: fundamental parameters — large-scale structure of 
universe — methods: statistical — surveys 



1. INTRODUCTION 

Determining the relationship between the spatial dis- 
tributions of galaxies and dark matter remains one of the 
central problems of theoretical and observational cosmol- 
ogy. The radial distribution of galaxies within their host 
dark matter halos dictates the correlation function on 
scales smaller than the virial radii of the largest dark 
matter halos (< l/i _1 Mpc). Therefore, the measured 
galaxy correlation function itself is a powerful tool for 
shedding light on how galaxies trace the underlying dark 
matter within halos. 

The halo occupation distribution (HOD) frame- 
work has been established as a robust method for 
modeling galaxy clustering by simply characterizing 
the biased relationship between g alaxies and mass 
(see, e.g.. iPeacock fc Smith! 120001; IScoccimarro et al.l 
200lt iBerlind fc Weinberg 2002; ICoorav fc Shethl [20021 : 
Zheng et al.ll2005[ ). The HOD describes this relation by 



has been used successfully to model gala xy clustering: 
on scales larger than r ~ lOO/t^kpc ( e.g.. TZehavi et al.1 



specifying the probability distribution P(N\M) that a 
halo of mass M contains N galaxies, along with a pre- 
scription for the spatial distribution of galaxies within 
halos. For this latter component, it is typicaly as- 
sumed that a single 'central' galaxy lives at the center 
of each halo, with additional 'satellite' galaxies tracing 
the dominant dark matter component. This assumption 
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2001 12005bl 201 It [Zheng et al.l [20071 [20091k " However, 
Watson et al.l (|2010l ) showed that this assumption does 
not work to explain the smaller-scale clustering of lu- 
minous red ga laxies (LRGs) in the Sloan Digital Sky 
Survey (SDSS; lYork et al.l [20001. In order to achieve a 
good fit to the clustering of L RGs on scales smalle r than 
~ 100/i _1 kpc (measured by iMasiedi et al.l 120061 ). they 
found that the radial profile of satellite LRGs must have 
a much steeper inne r slope compared to the Navarro- 
Frenk- WJrite_(NFWjNi^ dark matter 
profile. IWatson et al.l (|2010ft concluded that the distri- 
bution of satellite LRGs within halos is better described 
by an isoth ermal distribution. 

Recently, I Jiang et all ([201 ll ) measured the projected 
two-point correlation function w p (r p ) for several lumi- 
nosity sample s from the SPSS, e xtending the measure- 
ments done bv lZehavi et al.l (|2011l) on intermediate scales 
(~ 0.1 — 40/i~ 1 Mpc) down to extremely small galaxy- 
galaxy separations (0.01 — 7/i _1 Mpc). Our motivation 
for this paper is to model the innermos t data points 
0.01 — 0.4/i _1 Mpc) in the same vein as IWatson et al.l 
2010l hereafter W10) to see if there is a luminosity de- 
pendence of the radial profile of satellite galaxies. 

The pape r is lai d out as follows. In iJ5J we review the 
I Jiang et al.l (|2011l hereafter Jll) measurements. In £j3j 
we discuss our modeling method in the following manner: 
we provide an overview of the general technique 



used to model the small-scale correlation function; in 
we revisit the W10 four-parameter PNM model that al- 
lows the the probability distribution P(N\M) to vary, 
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and the five-parameter PNMCG model that also allows 
the radial distribution of galaxies within halos to vary. 
We present our results in 21 Finally, we summarize and 
discuss the implications of our results in fj5] Throughout 
the paper, we assume a ACDM csomology with Q m = 
0.25, il A = 0.75, n b = 0.04, h = 0.7, <j 8 = 0.8, n s = 1.0. 

2. DATA 

J 11 measured the projected two-point correlation func- 
tion down to very small scales (0.01 < r < 7/i~ 1 Mpc) 
for a large range of volume-limited galaxy luminosity 
threshold samples. Galaxy number denisities and me- 
dian redshifts for each sample can be found in Table [TJ 
These samples are constructed from the NYU VAGC 
(jBlanton et alJl2005[ ) V7.2 data (jAbazaiian et all 12009( 1 
which contains 8.6 x 10 5 SDSS Main Sample galaxies 
(|Strauss et al.|[2002h . Spectra for each luminosity thresh- 
old sample are cross-correlated with the full imaging sam- 
ple to compute w p (r p ) (details of the method can be seen 
below). The full imaging sample consists of ~ 10 8 galax- 
ies drawn from the SDSS imaging catalog. 

Measuring w p {r p ) on very small scales is a non-trivial 
task. The SDSS spectroscopic sample suffers from incom- 
pleteness due to the physical size of the fiber-optic cables, 
which impede the ability to take spectra of two galaxies 
closer than 55" on the sky. These "fiber collisions" thus 
result in a minimum pair separation of 55" . This effect is 
parti ally mitigated by the tiling method (jBlanton et al.l 
2003), which overlaps spectroscopic plates in order to 
yield full sky coverage. In plate overlap regions, galaxy 
pairs closer than 55" can be recovered. However, this still 
leaves ~ 9% of galaxies without measured redshifts (Jll) 
and this can strongly affect clustering measurements on 
very small scales. Followi ng the approach described in 
iMasiedi et all (|2006l . I2008D . Jll used a cross-correlation 
technique between the imaging and spectroscopic sam- 
ples to correct for fiber collisions and obtain unbiased 
w p (r p ) measurements down to the 10/i _1 kpc scale. De- 
tails of the method can be found in [j3]of Jll. 

A second potential problem that affects very small 
scales is that galaxies may be overlapping and the light 
within this region needs to be properly distributed be- 
tween the two galaxies. This is known as "deblend- 
ing" . M06 found a systematic error in the SDSS pipeline 
wherein too much light was allocated to the dimmer of 
the two galaxies. As a result, a galaxy that may have 
been too dim to make the LRG brightness cut could now 
be included in the sample. This increased the number 
of small-scale pairs and boosted the correlation function 
on very small scales. M06 quantified this effect and cor- 
rected their measurements accordingly. Since the physi- 
cal sizes of galaxies decrease rapidly with decreasing lu- 
minosity, this photometric deblending error diminishes 
in lower luminosity samples. For this reason, Jll ig- 
nored this effect in the luminosity samples they consid- 
ered (M r < —18 through —21). Nevertheless, we esti- 
mate the maximum effect of deblending by applying the 
M06 LRG correction to the M r < —21 sample and re- 
peating our analysis. We find that our results do not 
change qualitatively. 

We restrict ourselves to modeling the Jll data points 
from - 10 - 400/i -1 kpc. These can be seen in Fig- 
ure [1] which shows w p (r p ) scaled by an r~ x power law, 
with each panel corresponding to a distinct luminosity 



threshold. The bottom right panel of Figure [T] shows the 
LRG data points measured by M06 and modeled by W10. 
Along with these 9 data points used in our modeling, we 
also incorporate the measured number density for each 
luminosity sample (see Table [1} , providing a tenth data 
point. We use the full covariance matrices from the 9 
Jll w p (r p ) data points for each sample, which were es- 
timated using jackknifc resampling of 50 regions on the 
sky. The same jackknife samples were used to estimate 
the error on each calculated number density. 

3. REVIEW OF THE METHOD 

3.1. The HOD and the Galaxy 2PCF 

Here we briefly review the method used in W10, which 
is based on the halo occupation distribution (HOD) for- 
malism. The HOD fully characterizes the number, veloc- 
ity and spatial distribution of galaxies within dark matter 
halos. The probability distribution that a virialized dark 
matter halo of mass M will host N galaxies is designated 
as P(N\M). (N) M is the first moment of P(N\M), and 
is the mean number of galaxies as a function o f halo 
mass. Motivated by theory (e.g.. iBerlind et al.l 120031 : 



iKravtsov et alll2004l;izheng et alll2005D . we consider cen- 
tral galaxies that live at the center of their host halo 
and satellite galaxies that orbit within the host poten- 
tial as separate terms. We thus write the first moment 
of the HOD as (N)m = 1 + (N sat )M. Furthermore, we 
assume that there is some minimum halo mass, M m ; n , 
below which a halo will contain no galaxies and above 
which there will always be at least one central galaxy. For 
the satellite component, we adopt the parametric form, 

(N sa . t ) M = exp[-AV(A4alo - A/min)] X (Afhalo/Ml) a , 

where the satellite galaxies obey a power-law function of 
slope a with an exponential cut-off at the low mass end 
at Mo. Mi is the characteristic mass scale where a halo 
will contain, on average, one central and one satellite 
galaxy. 

To calculate the mean number of satellite galaxy 
pairs (N sat (N sat — l)),the second moment of the satel- 
lite P(N\M), we assume that the number of satellite 
galaxies in a halo of mass M follows a Poisson dis- 
tribution of mean (N sat )M- This sort of HOD pa- 
rameterizatio n is widely used t o model galaxy clus- 



terin g data (IZehavi et all l2005albl 12011 
20051: IConroy et all 12001 IZheng et al.l 
Watson et al.ll2011( l. 



Our fiducial model for characterizing the spatial dis- 
tribution of galaxies within their host halos places 
one galaxy at the center and assuming that the satel- 
lites trace the underlying dark matter density distribu- 
tion. The dark matter density profiles are described 
by the NFW relatio n p(r) cx (c^—)- 1 ^ + c^—)- 2 
(jNavarro et al.llT997j ). where c is the concentration pa- 
rameter, c = R v i T /r s , and r s is the characteristic scale 
radius. We use the virial definition of a halo to calcu- 
late the virial radius of host halos, such that i? v ; r = 
((3A/)/(47rA v i r p)) 1 / 3 , where A vir is the mean halo over- 
density (A v ir = 200), and p is the mean density of the 
universe. The con centration - mass relati on is given by 
I Zheng et all (|2007l ) for the modification of iBullock et all 
(f200U) 



: c 



_£q_ 



__ x (&g*.)-P, where c = 11, M, is the 
non-linear collapse mass at the median redshift of the 
sample for our choice of cosmology (M* is redshift de- 
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TABLE 1 

INPUT DATA FOR ALL LUMINOSITY SAMPLES 



M r 


n g 


^med 


log(A/„) 


-18.0 


3.209 (0.169) 


0.032 


12.4228 


-18.5 


2.405 (0.131) 


0.039 


12.4147 


-19.0 


1.689 (0.092) 


0.046 


12.4057 


-19.5 


1.326 (0.058) 


0.059 


12.3887 


-20.0 


0.749 (0.018) 


0.078 


12.3647 


-20.5 


0.377 (0.009) 


0.097 


12.3397 


-21.0 


0.123 (0.002) 


0.116 


12.3057 



Note. — n g is measured in units of 10 — 2 /i 3 Mpc~ 3 and the 
associated error from 50 jackknife samples on the sky is shown 
in parentheses. z med is the median redshift as measured in 
Jll. M* is the non-linear collapse mass and is in units of 
h~ 1 M Q . 



pendent and is thus uniquely denned for each luminosity 
sample and is given in Table [T]) and (3 = 0.13. 

The mean number density of galaxies can be calcu- 
lated for a given HOD by weighting the abundance of 
halos by (N)m and integrati ng over all halo mas s (see 
Eq.[2] of W10). We adopt the IWarren et al.1 l|2006f ) halo 
mass function, dn/dM, in all of our calculations, but our 
results are not sensitive to the specific choice of mass 
function. 

To model the galaxy two-point correlation function 
(2PCF), £ gg , we use the halo model. £ gg can be de- 
composed into contributions due to galaxies residing in 
the same halo (the "one-halo" term) as well as galaxies 
living in separate, distinct halos (the "two- halo" term). 
Therefore, £ gg can be writte n as the sum of the o ne- 
halo and two-halo terms (e.g.. ICoorav fc S hcth 20 03 fo r 
this particular form of the equation see IZhengl 12004): 

£gg( r ) = £gg al ° + £gg al ° + !• Since we are probing such 
small scales we need only consider the one-halo term for 
modeling £. However, it is possible that omitting the 
two-halo term entirely could cause a bias in the best- 
fit parameters for the one-halo term. This bias would 
be largest for the dimmest galaxy sample, because these 
galaxies live in the smallest halos, and the two-halo term 
will become important at small scales. However, even for 
these galaxies the effect is limited. For instance, while 
it is true that isolated —18 galaxies will reside in halos 
with radii much smaller than Q.4/i Mpc, the relevant 
halo size is set by M\. Furthermore, this is a threshold 
sample, so although the —18 sample will be dominated 
by galaxies of that brightness, there will still be much 
brighter galaxies that live in larger halos that influence 
Mi. As seen in the upper-left panel of Fig. [31 (Mi) for 
the —18 sample is 10 12 77 . This corresponds to a virial 
radius of order 0.4/j _1 Mpc. To carefully test the influ- 
ence of the two-halo term, we have have measured where 
the 1- and 2-halo terms cross for the —18 sample. This 
occurs at roughly l/i -1 Mpc, so the influence of the 2- 
halo term at 0.4/i _1 Mpc has a small effect on w p (r p ). 
We find that the amplitude of the 2-halo term is ~ 7% 
of the one-halo term at 0.4/i _1 Mpc for the —18 sample. 
Nevertheless, we do approximately account for the effect 
of the 2-halo term on our modeling, and we discuss these 
assumptions in detail below. 

The one-halo term depends on the second moment of 
P(N\M), as well as the pair separation distributions of 
central-satellite and satellite-satellite pairs (see Eq.[4] of 
W10 for the particular form). The pair separation dis- 



tribution for central-satellite pairs is essentially the same 
as the density profile of satellite galaxies. The satellite- 
satellite pair distribution is the convolution o f the density 
profile with itself and, as in W10, we use the lSheth et al.l 
(|2001[ ) calculation for the convolution of a truncated 
NFW profile with itself. Therefore, we can start to see 
the link between the shape of £(r) and the shape of the 
satellite profile on small scales. This is discussed in detail 
in SI 

For each luminosity threshold, Jll measured the pro- 
jected correlation function, so we convert our theoretical 
real-space correlation function £(r) to w p {r p ) (Eq.[5] of 
W10). At each luminosity threshold, w e integ rate up 
to the 7r max value given in IZehavi et al.l (|201lD . which 
ranges from 40 — 60/i _1 Mpc. The method of Jll results 
in a 7r max that is effectively larger than our choice. How- 
ever, the difference is insignificant at the small scales we 
model. Integrating out to large scales means that the 
two-halo term cannot be entirely ignored. Since calcu- 
lating the two-halo term correctly is fairly complex and 
its contribution to w p (r p ) is minimal at the scales of the 
data points that we model, we use a simple approxima- 
tion inst ead. We use a two- h alo term whose shape is the 
same as IZehavi et al.l ()2011l ). but whose amplitude can 
change with HOD parameters. Thus, for a given set of 
HOD parameters, we first calculate the large-scale bias 
of galaxies, b g , 

b 9 = n e 1 J™ dM^b h (M)(N) M , (1) 

where bh(M ) is the large-scale bias of halos from 
I Tinker et al.l ((2008) . The two- halo term has a simple re- 
lation to the matter-matter correlation function on large 
scales: £ gg al ° = b 2 g ^ mm . Thus, the relative two- halo terms 
of two different galaxy samples really only depends on the 
ratio of the bias of e ach sample. We cons ider the best-fit 
two-halo term from IZehavi et al.l ()2011l ) for each lumi- 
nosity and append a newly defined two-halo term to our 
model one-halo term, 

e2halo I ^9 \ 2 c 2halo / 

£gg = [bj ) %g ' ( 2 ) 

where primes designate the IZehavi et al.1 (|2011h values. 
In effect, we allow the amplitude of the 2-halo term to 
vary with HOD parameters, but we keep its shape fixed. 
The error introduced by this approximation is completely 
negligible as we have checked that even a shift as large as 
20% in the amplitude of the two-halo term has no appre- 
ciable effect on w p (r p ) at the maximum pair separation 
that we consider (0.4/i _1 Mpc) for any of the luminosity 
threshold samples. This test also emphasizes our insensi- 
tivity to the difference bet ween our choi c e of in tegrating 
out to the 7r max given by IZehavi et al.l (|2011l) and the 
technique used by Jll. 

3.2. The PNM and PNMCG Models 

W10 considered 4 distinct models for modeling the 
small-scale w p {r p ) of LRGs. In this paper we need only 
consider two of the models, defined as: 

• PNM - in this model the only free parameters are 
those associated with the P(N\M) distribution as 
described in g3TT] - M min ,M , M 1} and a. While 
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Fig. 1. — Model fits to the projected correlation function for luminosity thresholds spanning M r < —18 to LRGs. The points in the 
M r < —18 through M r < —21 panels show the Jiang et al. (2011) w v (r v ) measu rements (multiplied by r p ) and their associated jackknife 
errors. The bottom right panel shows results for LRGs from Watson ct al. (2010). The black curves show the best-fit PNM models, which 
use four-free parameters describing the probability distribution P(N\M) (M m j n , Mi, Mo, and a). The PNMCG model, shown in cyan, 
considers two additional free parameters (but Mo is fixed to M m ] n ): (1) / ga i, which relates the concentration of the density profile of 
satellite galaxies relative to dark matter (/ ga i = c ga i/c), and (2) 7, which allows for the inner slope of the density profile to differ from the 
dark matter distribution. The reduced x 2 values for each best-fit model are listed in each panel. 



these parameters are free to take on any value, the 
spatial distribution of satellite galaxies within ha- 
los is fixed to an NFW dark matter density pro- 
file. Since there are 10 data points (9 from the 
Jll w p (r p ) measurements along with the measured 
number density for each luminosity sample) and 4 
free parameters, the PNM model has 6 degrees of 
freedom. 



• PNMCG - this model allows for the same free pa- 
rameters as the PNM model, with the exception of 
Mq. Mq was unconstrained in our intitial MCMC 
runs (this was also the case in W10), thus we fixed 
Mq to M m ; n throughout our analysis. The PNMCG 
model also considers a parametrized density profile 
for satellite galaxies. As in W10, we allow both 
the concentration and the inner slope of the density 
profile to be free. The satellite galaxy concentra- 
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M r <-21.0 
M r < -20.5 




M r <-18.0 



0.01 0.10 
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Fig. 2. — The real-space correlation function residuals from 
an r~ 2 power law from the PNMCG best-fit models before con- 
verting to the projected correlation functions of Fig.[TJ Amplitudes 
have been arbitrarily shifted for clarity. The slope of £(r) on small 
scales is a reflection of the central-satellite pair distribution, which 
is essentially just the density profile itself. There is a strong lumi- 
nosity dependence of the slope of £(r) on small scales, becoming 
steeper and steeper for the brighter galaxy samples and this carries 
directly over to the luminosity dependence of the radial profile of 
satellite galaxies. The dotted lines show r^ 1 and r~ 2 power-law 
slopes. 

tion can differ from the dark matter concentration 
(denned in §3.ip through the free parameter f ga \, 

Cgal = /gal X C. (3) 

The inner slope of the density profile is no longer 
fixed to -1, as is the case for an NFW profile, but 
rather is specfied by the free parameter 7. As in 
W10, we adopt the following density profile for 
satellite galaxies 

p(r) = — (4) 

For an NFW profile 7 = 1, however, the PNMCG 
model allows 7 to take on any value from 0—4. This 
model considers the same number of data points as 
the PNM model, but now there are 5 free parame- 
ters (M m i„, Mi, a, /g a i and 7) resulting in 5 degrees 
of freedom. 

As detailed in §3.4 of W10, we use a Markov Chain 
Monte Carlo (MCMC) method to probe the parameter 
space for a given set of parameters (see iDunklev et al.l 
120051 for details on MCMC techniques). When a chain 
has converged, we can find the most likely value for each 
parameter by calculating the mean of its distribution. 
Errors for each parameter are given by the extrema of the 
middle 68.3% of the distribution. Best-fit parameters are 
found by the combination of parameter values for which 
X 2 is a minimum. 

4. RESULTS 



J 11 measured w p (r p ) to very small scales over a large 
range in luminosity thresholds. Their measurements 
nicely overlap those of lZehavi et al.l (|2011l ) on intermedi- 
ate scales (~ 0.2 - 7/i _1 Mpc, see Fig. 14 in Jll), and 
extend down to very small scales with the innermost 
data point for each sample at ~ 10/i~ 1 kpc. As discussed 
in 33. 2[ for each luminosity threshold sample, we model 
w p (r p ) with, (1) the PNM model, which only varies pa- 
rameters that determine the number of galaxies in a given 
halo, but forces the satellite galaxies to have an NFW 
spatial distribution within their halo, and (2) the PN- 
MCG model, which also allows the spatial distribution 
of satellite galaxies to vary within halos. 

Figure 1 shows our modeling results for each luminos- 
ity sample. w p (r p ) has been scaled by an r" 1 power law 
to more clearly highlight any discrepancies between the 
PNM and PNMCG models. Each panel shows the SDSS 
data points as well as the best-fit model for the PNM 
(black curve) and PNMCG (cyan curve) cases. It is clear 
from the figure that as we go to higher luminosities, the 
PNMCG model provides a significantly better fit to the 
data. We find that the P(N\M) parameter distributions 
are nearly the same for the two models, differing by, at 
most, ~ 3<T. Therefore, the improved fits for the PN- 
MCG model principally arise from the freedom to vary 
the density profile of satellite galaxies. Varying the den- 
sity profile is thus necessary to find a better fit to the 
data as we go to higher luminosities. 

We note that the reduced \ 2 values (listed in each 
panel) are in many cases quite high, even in the PN- 
MCG case. This could mean that the PNMCG model 
contains incorrect assumptions or does not have enough 
freedom. On the other hand, it could mean that the 
Jll jackknife errors are underestimated. To check the 
impact of the error estimates on our modeling, we re- 
estimated errors for the M r < — 20 sample using mock 
galaxy catalogs from the LasDamas project (McBride et 
al. in prep.). We used 160 catalogs 1 and measured the 
dispersion of w p (r p ) between the catalogs, using the same 
binning method as Jll. We then applied the fractional 
error (with respect to the mean of all mock measure- 
ments) to the data (non-mock) measurement to estimate 
the absolute errors and full covariance. Finally, using our 
new mock based error estimates, we re-ran the MCMC 
chains for the PNMCG model. We then compared the 
best-fit parameter values for the two fits and found that 
the parameters did not change significantly. By this, we 
mean the difference in \ 2 between the two best-fit points 
were within la when evaluated with either of the likeli- 
hood surfaces (from each of the two MCMC runs with 
different errors). We conclude two things from this test: 
(1) our somewhat high \ 2 values are not overly concern- 
ing, and are likely due to a slight underestimate of the 
errors from jackknife re-sampling on the data, and (2) 
this issue does not seem to affect any of our conclusions. 

We now investigate the luminosity dependence of the 
radial distribution of satellite galaxies and the degree 
to which it differs from an NFW distribution. As dis- 
cussed in 33.11 when constructing the real-space correla- 
tion function in the halo model, the one-halo term con- 
siders contributions from central-satellite and satellite- 

1 North-only SDSS footprint from the LasDamas "gamma" data 
release. 
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Fig. 3. — The radial profile of satellite galaxies as a function of luminosity. Each panel shows results for a specific luminosity sample, and 
for a halo mass that is chosen to be the mean value of M\ from the PNMCG MCMC chain (listed in the top right of the panel). The solid, 
yellow bands show the NFW profile, with a 20% uncertainty to account for the possible ambiguity in the dark matter distribution. Green 
curves show the profile corresponding to the best-fit PNMCG model. Each panel also contains 50 grey curves representing 50 randomly 
drawn links from the top 95% of the PNMCG MCMC chain, after sorting in \ 2 from lowest to highest values. These 50 curves thus span 
the 95% confidence region allowed by the data. As we probe higher luminosities, the radial profile of satellite galaxies strongly deviates 
from an NFW distribution on small scales, showing that luminous galaxies are poor tracers of the underlying dark matter. 
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Fig. 4. — Top Panel: The slope of the inner density profile of 
satellite galaxies as a function of galaxy luminosity. The filled red 
squares (and connecting lines) show the mean value of 7 from the 
MCMC chain of the PNMCG model for each luminosity sample 
and the green squares represent the best fit. The light blue and 
dark blue bands are the associated errors from the extrema of the 
middle 68.3% and 95% of the distribution, respectively. 7 = 1 
corresponds to an inner slope for an NFW dark matter profile, and 
is shown as a solid black line. The dotted black lines highlight an 
assumed 20% inaccuracy in the dark matter profile. Bottom Panel: 
The same procedure as the top panel, but for the parameter / ga i 
which relates the galaxy and dark matter concentrations (/ ga i = 
c gal/c). /gal and 7 are intrinsically linked and there is a strong 
trend towards inner slopes becoming steeper than NFW as we go 
to higher luminosities, with corresponding decreasing values of / ga i- 



satellite pairs. The central-satellite contribution, which 
is essentially just the density profile itself (see Eq. 4 of 
W10), is steeper than the satellite-satellite pair contri- 
bution and thus dominates the correlation function on 
the very small scales that w e are considering (e.g., Fig- 
ure 4. of lZheng et al.1 12(309"). Therefore, the luminosity 
dependence of the slope o/£(r) on small scales can give a 
direct indication of the luminosity dependence of the ra- 
dial profile of satellite galaxies (though the slope of 
will be less steep than the slope of the radial density pro- 
file due to the dampening effect of the satellite-satellite 



term). Figure [2] shows the residuals from an r~ 2 power 
law from the PNMCG best-fit models before converting 
to the projected correlation functions shown in Figure [TJ 
The amplitudes of the curves have been arbitrarily stag- 
gered simply to make the plot more clear. The dotted 
lines highlight the cases of r~ l and r~ 2 power laws. The 
slope of £(r) is clearly a strong function of luminosity, 
being close to -1 for low luminosity galaxies and going 
more and more towards -2 for the M r < —21 sample and 
even steeper for LRGs. The W10 result for the steepness 
of the sl ope of £(r) for LRGs on small sc ales was also 
found by lAlmeida et all (|2008[ ) . Using the iBower et all 
(2006) s emi-analytic model a pplied to the Millenium sim- 
ulation (|Springel et al.l [20051) . they found that the LRG 
real-space correlation function follows an ~ r ~ 2 -07 p 0wer 
law shape down to the ~ 10/i _1 kpc scale. 

We next wish to directly investigate the radial profiles 
of satellite galaxies that are required by the data and 
compare them to the NFW profile. For each luminosity 
sample, we choose a halo mass equal to the mean value 
of Mi in the PNMCG MCMC chain for that sample. We 
choose Mi because it represents the typical size halo that 
contributes central-satellite pairs to £(r) (smaller halos 
have no satellites and larger halos are rare). Specifying 
the halo mass sets the amplitude, virial radius, and dark 
matter concentration in Equation 21 The radial profile 
then only depends on the parameters / ga i and 7. We then 
take the full MCMC chain for the PNMCG model and 
sort it by x 2 from lowest to highest. We then randomly 
draw 50 links from the top 95% of the chain. Each of 
these links has distinct values of / ga i and 7, which we 
insert into Equation [4] in order to construct individual 
radial profiles. The light grey curves in each panel of 
Figure [3] show these 50 profiles (multiplied by 4nr 2 dr 
to convert them from density into mass profiles). These 
curves thus span the 95% confidence region allowed by 
the data. The green curve in each panel shows the radial 
profile corresponding to the best-fit PNMCG model. To 
compare with the NFW profile for dark matter, we also 
plot the case / ga i = 1,7 = 1, shown by the yellow bands. 
We assume that the dark matter profile of halos has a 
20% uncertainty, which we represent as the thickness of 
the bands (see Sj5] for discussion of the uncertainties in 
the dark matter distribution in halos). For reference, the 
scale radius for any particular sample occurs where the 
dark matter radial profile "turns over". Figure [3] shows, 
once again, that galaxies of increasing luminosity deviate 
more from an NFW profile. In fact, the figure seems 
to suggest that there is a transition point, somewhere 
between an absolute r-band magnitude of -19.5 and -20, 
where the radial profile of satellite galaxies goes from 
being consistent with NFW to completely inconsistent. 
For the highest luminosity samples, the profiles approach 
a power-law shape. 

Figure 2] shows the luminosity dependence of the inner 
slope and concentration of satellite galaxies as found by 
the PNMCG model. The dark blue and light blue bands 
in each panel represent the middle 95% and 68.3% of 
the MCMC chains after first sorting in 7 (top panel) 
and sorting in / ga i (bottom panel). The filled red 
squares show the mean of the respective parameters in 
the MCMC chain for each luminosity sample and filled 
green squares show the best-fit values. The solid black 
line highlights the NFW dark matter inner slope of the 
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Fig. 5. — The 7 — / ga i parameter space as a function of luminosity, highlighting the relationship between the inner slope and concentration 
of the density profile of satellite galaxies relative to dark matter. For each panel, the dark blue and light blue regions are defined as the 
top 95% and 68.3% of the MCMC chain for a given luminosity sample after sorting in \ 2 ■ There are two points of reference designated as 
filled circles - the 7 = 1, / ga i = 1 combination that represents an NFW distribution, and the 7 = 2, / ga i — > combination that represents 
a singular isothermal sphere distribution (SIS). Red squares show the mean values of 7 and / ga i from our MCMC runs and green squares 
show the values for the best-fit model. The data for M r < —20 and brighter galaxies strongly rule out the NFW model and are in good 
agreement with an SIS distribution for M r < —21 and LRGs. 
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Fig. 6. — A more physically useful definition of concentration: 
Mi/io, denned as the fraction of satellite galaxies (or mass) that 
are enclosed within one tenth of the virial radius of a halo (see 
Eq. [5}. Top Panel: For each luminosity sample, we choose a mass 
equal to the mean value of Ah from the PNMCG MCMC chain, 
and we compute -Mi/iq for each link in the chain. Red squares 
(and connecting lines) show the mean of the Mj/kj distribution for 
a given luminosity, and errorbars for M 1/ / 10 show the extrema of 
the middle 68.3% of the distribution. Black points show M 1/ / 10 for 
an NFW distribution. Bottom Panel: Ratio of Afi/io for satellites 
with respect to dark matter. The dotted line highlights a ratio 
of unity. Satellite galaxies clearly become much more centrally 
concentrated than dark matter with increasing luminosity. 

density profile which corresponds to 7 = / ga i ~ 1 and 
the dotted lines represent the possible 20% ambiguity of 
an assumed NFW inner slope and concentration. The 
M r < —18, —18.5, —19 and — 19.5 samples have large 
spreads in 7 and / ga i values, as there are many parame- 
ter combinations that can yield a similar goodness-of-fit 
to the data. For M r < —20 and brighter galaxies, both 
7 and /gai have tighter distributions and deviate from 
NFW. Figure [4] shows that there is a strong luminosity 
dependence of the satellite galaxy profile, i.e. 7 is an in- 
creasing function of luminosity, and / ga i is a decreasing 
one. However, notice the unique coupling of / ga i and 7 
for the —20.5 case, where / ga i seems to deviate from the 
aforementioned trend. 7 for —20.5 satellite galaxies also 
strays from the trend, but it is still significantly larger 
than NFW, and still lies within la of the M r < -20 and 
M r < —21 values (/ ga i is within 2a). However, the rea- 
son for this outlier is unclear, though it may be due to 
the fact that w p (r p ) has a strong feature at roughly the 



scale radius of the typical host halo (see the M r < —20.5 
panel of Fig. [1]). This may cause a unique interplay 
between / ga i and 7, resulting in a mild departure from 
the trends. For M r < —20 and brighter galaxies, the 
mean value of 7 ranges between 1.6 - 2.1, significantly 
steeper than NFW. This range of 7 values persists all 
the way down to M r < —19 when considering the best- 
fit values. As mentioned in £[2] M06 found a systematic 
error in the manner that light in the overlapping region 
of LRG pairs was being allocated in the SDSS pipeline. 
This caused a slight boost in the small-scale correlation 
function and was corrected for by M06. This effect is 
expected to rapidly diminish for lower luminosity sam- 
ples, because the intrinsic sizes of galaxies will be smaller. 
However, we apply the same LRG correction to the Jll 
M r < —21 sample to model the maximum effect that 
this could have. While we expect (7) to decrease as a 
result of the innermost data points being scaled to lower 
w p (r p ) values, we still find that (7) ~ 1.6. Thus, even af- 
ter assuming a drastic over-correction, 7 is still strongly 
discrepant from NFW. 

Figure 0] showed the luminosity trends for 7 and / ga i 
by sorting the MCMC chains separately for each param- 
eter. However, 7 and / ga i are intrinsically coupled to 
one another so it is important to investigate the joint 
7 — /gai parameter space. Figure [5] Shows the 7 — / ga i 
parameter space at each luminosity, with the dark blue 
and light blue regions defined as the top 95% and 68.3% 
of the MCMC chain, after sorting in x 2 . There are 
two points of reference designated as filled circles: the 
7 = 1, /gai = 1 combination for an NFW distribution 
and a 7 = 2, / ga i — > combination representing a sin- 
gular isothermal sphere (SIS) distribution. Red squares 
represent the mean values of 7 and / ga i from our MCMC 
runs and green squares show the best-fit combination. As 
expected, the M r < —18 through —19.5 exhibit a broad 
range in 7 — / ga i combinations, with the NFW combina- 
tion lying within the top 68.3% region (with the excep- 
tion of the -18.5 sample). The 7 — / ga i regions clearly 
drift away from NFW for M r < — 20 and greater lumi- 
nosities and are well described by an SIS distribution for 
both the —21 and LRG samples. Figure [5] shows that the 
inner slope and concentration parameters are strongly 
degenerate with each other, particularly at low luminosi- 
ties. This can be caused by having the scale radius (the 
radius within the halo where the slope transitions from 
—7 to -3) too close to the innermost data point. When 
this happens, the data cannot accurately constrain the 
inner slope, and only constrains 7 and / ga i through their 
degenerate contribution to the amplitude of the density 
profile. A simple calculation of the scale radius using the 
virial definition of a halo (see ij3.1|) gives r s ~ 10/i _1 kpc 
for the M r < — 18 sample (adopting a halo mass equal 
to the mean value of Mi ) , which is precisely at the scale 
of the innermost data point of w p (r p ). This leads to the 
large spread in 7 — / ga i combinations. Of course, as we 
go to brighter samples, the scale radii increase, and the 
parameters are better constrained. 

For each luminosity, as the inner profile steepens, / ga i 
decreases (with the slight deviation for the M r < —20.5 
case). The decrease in concentration does not necessarily 
imply that the concentration of satellites is low in the 
traditional sense - i.e., that there are fewer galaxies at 
small radii - but rather that the radial profiles simply 
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evolve toward a more power-law form. For example, as 
/gal — > in the case of LRGs, this does not mean that the 
LRG satellites are less concentrated than the underlying 
dark matter. As / ga i — > 0, the scale radius is pushed 
far outside the virial radius, resulting in a pure power- 
law distribution of slope —7. Since the value of / ga i can 
be misleading when thinking about the concentration of 
satellite galaxies, we consider an alternative definition 
of "concentration" that is more physically useful: the 
fraction of satellite galaxies (or mass) that are enclosed 
within a tenth of the virial radius i? v ir- We call this 
quantity Mi/ W and compute it by integrating any radial 
profile out to one-tenth the virial radius of the host halo 
and dividing by the total mass out to i? v ir 

M 1/w = M } r <°- 1R ^\ (5) 
1/10 M(r < R vir ) K ' 

Figure H] (top panel) shows M 1 / w as a function of lu- 
minosity. For each luminosity sample, we assume a halo 
mass equal to the mean value of Mi from the MCMC 
chain, and we calculate M 1 ^ 10 for every link in the chain 
using each link's values for 7 and / ga i. We then find the 
mean of the distribution of M 1 i w values (denoted by the 
red squares and connecting lines). Errors for M 1 / w are 
given by the extrema of the middle 68.3% of the distribu- 
tion after sorting by Mi/i . The black filled circles show 
My 10 f° r a dark matter NFW profile. The dark matter 
concentration drops with luminosity because the mean 
value of Mi increases with luminosity a nd concentra- 
tion i s a decreasing function of halo mass (| Bullock et al.l 
I2001D . The bottom panel of Figure |5] shows the ratio 
of Mi/10 f° r satellites with respect to dark matter. The 
figure clearly shows that low luminosity satellite galax- 
ies (M r > —19.5) have spatial concentrations similar to 
dark matter, whereas the profiles of luminous satellite 
galaxies (M r < —20 and brighter) are ~ 2.5 — 4 times 
more concentrated than NFW. 

5. SUMMARY & DISCUSSION 

Do satellite galaxies trace the underlying dark mat- 
ter distribution? Our modeling of the small-scale w p (r p ) 
over an enormous range in galaxy luminosity thresholds, 
from SDSS Main galaxies to LRGs (M r < -18 through 
< —23), has revealed a strong luminosity dependence 
of the radial distribution of satellite galaxies. We have 
found that for the low luminosity samples (lower than 
M r < —20), there is a wide range of satellite galaxy con- 
centrations (equal to / ga i times the NFW dark matter 
concentration) and inner density profile slopes, —7, that 
are consistent with the data. This lack of constraining 
power is due to a strong degeneracy between / ga i and 7, 
which, in turn, is possibly due to the fact that the small- 
est scale data point in w p {r p ) is roughly around the scale 
radius for halos hosting galaxies of this size. Neverthe- 
less, low luminosity satellite galaxies have radial distri- 
butions that are generally consistent with the NFW dark 
matter distribution. When we move to higher luminos- 
ity galaxies (M r < —20 and brighter samples) however, 
satellite profiles change dramatically, with the most strik- 
ing feature being that 7 jumps to much higher values. 
Our results show that for M r < — 20 and brighter galax- 
ies, 7 ranges from ~ 1.6 — 2.1, highly discrepant from 
NFW, even after assuming a 20% inaccuracy in the dark 
matter profile. Luminous satellite galaxies are thus poor 



tracers of the underlying dark matter within halos on the 
small scales that we probe. Since the effects of / ga i and 7 
on the radial profile are intertwined, we have calculated 
a more physically useful quantity Mi/i , which gives the 
fractional amount of mass enclosed within one-tenth of 
the virial radius. We find that Mi/ 10 is a strong function 
of luminosity, being consistent with NFW for low lumi- 
nosity galaxies and being ~ 2.5 — 4 times more concen- 
trated than NFW for galaxies brighter than M r < —20. 

Several other studies have also investigated the radial 
distribution of galaxies within halos. When consider- 
ing faint satellit e s arou nd ~ L* galaxies, no consensus 
emerges. iChenl (|2009l ) foun d the satell i te dis tribution 
to be consistent with NFW. iMore et al.l ()2009t ) found a 
cored satellite profile best descibed by a 7 = 0, ,f ga i = 
0.5 co mbination, inconsistent with NFW. Carlbcr g et al.l 
(2009) found that the radial profiles of dwarf satellites 
around nearby galaxies are much more concentrated 
than NFW, rendering all of these studies in disaccord. 
Furthermore, the discrepancies persist when considering 
a bro ader range of galaxy systems. iNierenberg et all 
(|2011l) recently studied the radial profile of satellites 
around massive, early-type galaxies at intermediate red- 
shifts. Assuming a power-law profile model, they find 
that the s atellites h a ve an isothermal distribution of 
slope -1. IGuo et al.l (|2012t) considered SDSS satellite 
galaxies around host central galaxies of a wide luminos- 
ity range and, in general, satellite number density pro- 
files were shown to be consistent with NFW. The NFW 
dark matter concentration was found to decrease with in- 
creasing satellite galaxy luminosity, independent of the 
host galaxy luminosity (except for the brightest centrals) . 
They also detected a slightly steeper profile for fainter 
satellites. These varying results imply that uncovering 
the spatial distribution of satellite galaxies is a difficult 
problem. In the end, we have found that there is a strong 
luminosity dependence of the radial distribution of satel- 
lite galaxies, wherein ~ L* and brighter satellites within 
group and cluster sized halos have a substantially steeper 
radial profile than dark matter on scales smaller than 
~ 100/i _1 kpc. We emphasize that while the scales in- 
volved in all of these studies are similar to ours, galaxy 
sample selection is different in each case, making it diffi- 
cult to directly compare their results with ours. 

Our results have shown that luminous galaxies are poor 
tracers of NFW. We now address the possibility that 
NFW is a poor model for the underlying dark matter 
distribution. In the established ACDM concordance cos- 
mology, recent studies invoking high resolution N-body 
simulations have shown that the mass profiles of ACDM 
halos slightly, but systematically, d eviate from NFW, 
becoming shallow er at smaller radii ()Stadel et al.ll2009t 
iDel PopolofeOlOD . In fact, they may be better described 
by Einasto profiles, which include a parameter a that 
controls how the logarithmic slope will vary with radius 
to accurately account fo r the fact that ha l o profiles seem 
to not be self-simi l ar (|Gao et al.l 120081 : iNavarro et al.l 
[20101 iLudlow et all l2010h. These studies find, on av - 
erage. 7 < 1 (|Graham et al.l 120061; INavarro et alll2010l) . 
However, on average, the simpler, two-parameter NFW 
model, which has a character istic r~ x inner slope, is ac- 
curate to within 1 - 20 % (|Bensonl 120101) . Moreover, 
iMandelbaum et all (|2008l ) fit galaxy cluster weak lens- 
ing profiles on small scales and found that NFW and 
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Einasto profiles gave the same result to within several 
percent. Therefore, for simplicity, and to be consisten t 
with previous modeling of LRGs bv lWatson et all (2010), 
we assumed an NFW distribution throughout our model- 
ing analysis. We conclude from the aforementioned stud- 
ies that NFW is not a poor model for the dark matter 
distribution in collisionless simulations and, to the extent 
that it is, the true profile is even more discrepant from 
our results for luminous satellite galaxies. 

Of course, by assuming a pure dark matter profile es- 
tablished from high resolution N-body simulations, the 
effects that baryons can have on dark matter are not con- 
sider e d. Adiabatic Contracti o n (AC) (iBlumenthal et ah! 
fl98l iRvden fc Guml [l987t iGnedin et al.l 12004ft may 
cause a steepening of the inner slope of the dark 
matter density profile as the gas condenses and sinks 
to to the center of the dark matter potential well 
dDiemand et alj|2004l iFukushige et ai1l2004 i Reed et al.l 
l2005HDel Popolo fc Kroupall2009| ) . However, the major- 
ity of results indicate that the steepening is insubstantial 
or ca n vary widely from halo to halo (e.g.. [Tissera et all 
|2010() . In fact, while AC may cause a steepening of 
the inner profile, over time, major and minor merg- 
ers can cause the dark mat t er profile to become shal- 
lower (lEl-Zant et all l2001t iRomano-Di'az et al.l 120081 
l2009t I Johansson et al.1 l2009f ) . Observational confirma- 
tions of AC are notoriously difficu lt. In the case 
of gal axy clusters, X-ray analysis by IZappacosta et al.l 
(2006) snowed that processes during halo formation (e.g. 
gravitational heating from merger events) counteract AC 
and the mass profiles wer e still well described by NFW. 
iMandelbaum et all (|2006| ) used galaxy-galaxy lensing to 
show that the mass density pr ofile of LRG c l usters is con- 
sistent with NFW. Recently, ISchulz et all (|2010l ) stud- 
ied the profiles of a large sample of SDSS ellipticals and 
found that their dynamical mass measurements suggest 
that the measured excess mass on small scales may be 
consistent with the AC hypothesis. However, for bright 
galaxies, our results are still not reconcilable with cur- 
rent AC models. We thus conclude that luminous satellite 
galaxies are indeed poor tracers of the underlying dark 
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matter distribution, even accounting for the effects of 
baryons on the dark matter. 

It is perhaps not surprising that galaxies do not be- 
have like test particles as they orbit within a dark mat- 
ter potential well. Being massive and extended, they 
are subject to dynamical mechanisms that would not af- 
fect test particles, namely dynamical friction and mass 
loss due to the tidal field of their host halo. In fact, 
these mechanisms should affect large and massive ob- 
jects more than small ones, which could explain the lu- 
minosity trend that we observe. Satellite galaxies are 
thought to reside within subhalos, so it makes more sense 
to compare the radial distribution of galaxies to that of 
subhalos, rather than just dark matter. While subha- 
los tend to have a less concentrated radial distribution 
than dark matter dGhigna et al.lll998l 120001: jGao et al.1 
l2004HDe Lucia et al.ll2004HNagai fc Kravtsovll2005[) . the 
subhalo distribution has yet to be extended down to the 
extreme small scales that we have probed. High reso- 
lution N-body simulations are just now becoming avail- 
able that allow for subhalos to be distinguished at the 
10/i _1 kpc separation level (e.g.. iKlypin et all I2010T) for 
all of the luminosity samples that we have studied. Of 
course, subhalos may not be perfect tracers of galaxies, 
but it will be interesting to sec whether the distribution 
of subhalos as a function of mass follows the same trend 
that we have uncovered for satellite galaxies as a func- 
tion of luminosity. If the primary cause of this trend is 
dynamical in nature (e.g., dynamical friction), it should 
also show up with subhalos. Either way, such a compar- 
ison will yield insight into the relation between subhalos 
and galaxies on very small scales, as well as the compli- 
cated non-linear processes occuring towards the centers 
of host halos, crucial for galaxy formation theory. 
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